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Abstract 

The high-temperature expansions of the spin-spin correlation function of the two-dimensional 
classical XY (planar rotator) model are extended by two terms, from order 24 through order 26, 
in the case of the square lattice, and by five terms, from order 15 through order 20, in the case 
of the triangular lattice. The data are analyzed to improve the current estimates of the critical 
parameters of the models. We determine f3 c = 0.5599(7) and a = 0.499(5) in the case of the square 
lattice. For the triangular lattice case, we estimate (3 C = 0.3412(4) and a = 0.500(5). 
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We have recently extended^, from the 21th to the 24th order, the high-temperature 
(HT) expansions for the spin-spin correlation function of the classical two-dimensional XY 
Berezinskii-Kosterlitz-Thouless(BKT) lattice model^ with nearest-neighbor interactions on 
the square(sg) lattice and briefly updated the series analysis for the susceptibility and the 
second-moment correlation-length. Here we present a further extension of the HT expansions 
from the 24th order through the 26th order for the sq lattice and from the 15th order to the 
20th order in the case of the triangular (tr) lattice. In our derivation of the series coefficients, 
we have used a non-graphical recursive algorithm^ based on the Schwinger-Dyson equations 
for the spin-spin correlations of the XY-model. The necessary computations were performed 
using a few nodes of a pc cluster for an equivalent single-processor CPU-time of seven 
weeks in the case of the tr lattice and a CPU-time ten times as long in the case of the 
sq lattice. Arguing as in Ref.[l] that the "effective length" (namely the ability to provide 
accurate numerical information) of the HT expansions of a given system on different lattices, 
is proportional to the computer time used in the series generation, we expect that the 
" effective length" of our expansions for the tr lattice, is only slightly smaller than that of the 
sq lattice series and moreover that their convergence will be faster, because the tr lattice is 
closely-packed. 

Let us recall that the lattice XY spin model with nearest-neighbor interactions is described 
by the Hamiltonian 

H{v} = -2J #(r) ■ v(r') (1) 

<nn> 

where v(f) is a two-component classical unit vector at the site r and the sum extends to the 
nearest-neighbor sites of the lattice. 

We have calculated the spin-spin correlation function 

C(0,x;(3)=<v(0)-v(x)>, (2) 

as series expansion in the variable (3 = 2J/kT, with T the temperature and k the Boltzmann 
constant, for all values of x for which the HT expansion coefficients are non-trivial within 
the maximum order reached. In terms of these quantities, we can form the expansions of 
the l-th order spherical moments of the correlation function: 

m^{(3) = J2\x\ l <v{6)-v{x)> (3) 

X 

and, in particular, of the reduced ferromagnetic susceptibility x(P) — m^(/3). 

The second-moment correlation length is defined, as usual, in terms of m^(P) and x(P)'- 

e(P) = m {2) (P)/4 X (P). (4) 

After completing our computations, we have been kindly informed by H. Arisue 6 that, 
in the particular case of the sq lattice, he has recently obtained, using the "finite lattice" 
method, remarkably longer HT expansions: through the 38th order for the nearest-neighbor 
spin-spin correlation C(0, 0; 1, 0; 0) and through the 34th order for the susceptibility and the 
second and fourth moment of the correlation function. 

When also these further extended data will be published, the agreement, through their 
common extent, between both sets of results independently obtained by completely different 
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methods, will be a significant check of correctness in consideration of the notoriously intricate 
nature of the high-order HT computations. 

In order to determine the critical parameters of the models, the HT expansions of the 
above defined quantities should be confronted to the following main predictions^ of the 
BKT renormalization-group analysis of the XY system. 

As (3 — > (3 C , the divergence of the correlation length is expected to be dominated by the 
characteristic singularity 

m ~ Zas(P) = exp(6r- CT )[l + 0(r)} (5) 

where r — 1 — (3/ (3 C and the universal exponent a is expected to take the value a — 1/2, 
while b is a non-universal (and thus lattice dependent) positive constant. 
The critical behavior of the susceptibility is predicted to be 

X{P) ~ C 2 rW = exp ((2 - n)br-) [1 + 0(r)] (6) 

The parameter rj represents the correlation-function exponent and is predicted to take the 
value r] — 1/4. 

At the critical inverse temperature (3 = f3 c , the asymptotic behavior of the two-spin 
correlation function as \x\ = r — > oo is expected 8 to be 

< *3).i >~(!^[1 + 0(^)1 (7) 

The value 9 = 1/16 is predicted for the second universal exponent characterizing the 
critical correlation function. 

Following a recent renormalization-group analysis^ of the two-dimensional 0(2) nonlinear 
cx-model, which should belong to the same universality class as the XY model, we shall 
assume that the critical behavior of x{P) does not contain singular multiplicative corrections 
by powers of ln(£) (or equivalently singular multiplicative corrections by powers of r). The 
possible existence of such corrections has long been numerically investigated^^ 1 ^ with 
conflicting and essentially inconclusive results. As we will indicate below, our series analysis 
is completely consistent with the conclusions of Ref. [§]. 

Let us now come to our results for the sq lattice: we have added the two terms 



376988970189597090587^ 2 , 62337378385915430773643 n2f , 

y(B) = ...-{ /r 5 H /r 6 + ... (8) 

Ay 1 384296140800 31384184832000 v ; 

to the expansion of the susceptibility. 

To the expansion of the second moment of the correlation function, we have added the 
terms 

5412508223507386985733313^ 7139182711315236460182251 ^, 

vcP'iB) = ... H /r 5 H (3 26 + ... (9) 

yHJ 13450364928000 7846046208000 v ; 

The lower-order coefficients of the sg-lattice expansions of x(P) an d of m^(P) have been 
recently tabulated in Ref.[l[ and therefore they will not be reproduced here. 

On the contrary, we have listed in Table d also the coefficients already known^^ through 
15th order in the case of the tr lattice, in addition to the five coefficients that we have 
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recently calculated for the nearest-neighbor spin-spin correlation, the susceptibility and the 
second moment of the correlation function. 

Let us now update the analysis of Refs.[l|J3] for the sq lattice expansions by including 
all coefficients that we have so far derived. The expected form eq. ((SD of the critical 
singularity suggests that the critical parameters should be conveniently obtained by an 
inhomogeneous-differential-approximant^ (DA) study of the location and the exponent of 
the leading singularity of the quantity ln 2 (x). This conclusion is supported also by a simple 
comparison of the distribution of the singularities on the real /3-axis for the DAs of x{0) an d 
of ln 2 (x) which, independently of theoretical prejudice, suggests that the analytic structure 
of the latter form is more suitable to a study by DAs. Notice that the choice of this 
particular function of xW) does not imply any biasing of a, as long as we do not make any 
selection of the singularities of the DAs for ln 2 (x) on the basis of their exponents. In our 
previous note 1 , we had already observed that a possibly more accurate analysis might be 
based on the obvious remark that, near the critical point, by eq. flSJ), we have to expect 
ln(x) = c sq /r a + c 2 + ... Assuming a ~ 1/2, a simple fit to the asymptotic form of ln(x), 
can determine c 2 — —1.5. We are then led to study also the simple generalization of the 
function ln 2 (x), defined by L(a,x) = (a + ln(x)) 2 , where a is some constant. The relative 
strength of the t~° and r~ 2cr singularities in the function L(a, x) is determined by the value 
of a. Therefore it might be numerically convenient to analyze the function L(a, x) with 
a ~ — c s 2 q , rather than simply ln 2 (x), because the dominant singularity of L(— c%*,x) should 
then be r~ 2a , with a small or vanishing t~° contribution and thus approximately a simple 
pole, if ex — 1/2. We can therefore expect that analyzing the DAs of L(—c 2 q , x) wm make it 
possible to determine with good accuracy not only the position, (which is not very sensitive 
to the choice of a), but also the exponent of the critical singularity. 

In our analysis of the HT series, we have restricted to a class of quasi- diagonal second- 
order DAs quite similar to that used in Ref. [l|, namely to the [k, I, m; n] DAs defined by the 
conditions: 20 < k + l + m + n < 24 with k > 6; I > 6; m > 5 and such that \k— 1\, \l — m\ < 2, 
and 1 < n < 7. For the shorter tr lattice series, these conditions have to be modified in an 
obvious way. We have always made sure that our final estimates, within a small fraction of 
their spread, are independent of the precise definition of the DA class examined. 

In the case of the sq lattice, from an unbiased DA analysis of ln 2 (x) we obtain (3 C = 
0.5599(7), in complete agreement with the results of o ur p revious analysis of the 24-th order 
series and also with the MonteCarlo analysis of Ref. [1 1| which yielded (3 C = 0.55995(15). 



Notice that the uncertainty of this estimate was not given explicitly in Ref. [11], and that 



the value that we have nevertheless indicated is only a reasonable guess. We have obtained 
the estimate a = 0.499(5) for the singularity exponent by DAs of L(—c 2 q , x), chosen in the 



same class defined above and biased with the value of (3 C = 0.55995 given in Ref. [l l] . Its 



uncertainty accounts also for the spread in the estimate of the critical temperature used to 
bias the DA calculation. These results are illustrated in Fig{TJ 

In FigfSJ we have compared how the estimate of a depends on the value f3^ as used to 
bias the DAs, when either the quantity L(—c 2 q ,x) or l n2 (x) ^ s analyzed in a vicinity of 
pbias _ ^ _ Q.55995. Clearly the results obtained from the analysis of ln 2 (x) are much 
more sensitive to the bias value. 

Taking advantage of our extension of the HT series, we can also get some hint of the 
analytic structure of the susceptibility in a vicinity of the origin of the complex (3 plane. 
In Figj3l we have reported a scatter-plot showing the distribution of the nearby complex 
singularities for a large class of inhomogeneous first-order DAs of L(—c 2 q ,x)- We have 
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discarded only a few spurious real singularities with modulus less than 0.9/5 c . We can 
notice that, in addition to the physical singularities at ±/3 c and to many randomly scattered 
complex singularities, which we believe to represent only numerical noise, there is evidence 
of five pairs of complex-conjugate singularities, located just beyond the convergence disc of 
the series. They are indicated by clusters of DA singularities, which are likely to coalesce 
around tips of cuts with increasing series order. These results, which do not essentially 
depend on the value of the constant a in L(a,x), are reminiscent of the regular pattern of 
the nearby unphysical singularities (cuts) which were exactly determined for the ^lattice 
two-dimensional O(N) a-model 14 in the large N limit and conjectured^ to exist also for 
finite N>3. 

Let us now perform a similar update of the analysis for the tr lattice series. As shown in 
FigJU an unbiased DA study of ln 2 (x) yields a critical inverse-temperature (3 C = 0.3412(4), 
sizably improving the precision of our 14th-order estimate (3 C = 0.340(1) reported in Ref.[4|, 
but disagreeing with the DA estimate f3 c = 0.33986(4) obtained in the 15th-order study 
of Ref. [l2j. From an analysis of DAs of ln 2 (x) biased with the critical temperature (3 C = 
0.3412(4), we are led to the estimate a = 0.53(2). If, however, in analogy with the s^-lattice 
analysis, we consider the quantity L(a,x) with a = — c| r ~ 1.1, obtained from a fit of the 
asymptotic behavior of ln(x) as indicated above, we arrive at the estimate a = 0.500(4). 
It is also particularly interesting to notice that the pattern of the singularities of L(a, x) in 
the complex /3-plane is much simpler than that observed in the case of the sq lattice and 
consists of a single pair of complex-conjugate clusters, located farther from the border of the 
convergence disc than in the case of the sq lattice. This is an indication that the convergence 
of the tr-lattice HT series will be faster than in the sq lattice case, and so the argument that 
the effective lengths of the HT series for the two lattices are comparable, receives further 
support. 

In spite of the sizable extension of the series, both in the sq and in the tr lattice cases, 
the study of the indicator function H{(3) = ln(l + m( 2 )/x 2 )/l n (x) = 2^ + 0(r°") (or of anal- 
ogous functions of different correlation-moments), does not show very sharp improvements 
in the accuracy of the determination of the exponent 77, for which we consistently obtain the 
estimate 77 = 0.25(2). 

Assuming a = 1/2, essentially the same estimate of 77 is obtained from a DA analysis 
of the series expansions of r CT ln(x) = (2 — 77)6 + 0(r CT ) and of r cr ln(^ 2 //5)) = 2b + 0(r CT ) 
both in the case of the sq and of the tr lattices. Assuming 77 = 1/4, we can estimate the 
non-universal parameters h sq = 1.77(1) and bt r = 1.70(1). 

In order to get alternative estimates of rj, we may also try to analyze directly the large- 
order behavior of the expansion coefficients of x(P) an d of £ 2 (/3). Assuming a = 1/2 and 
using the known value of /3 C , we can fit the series coefficients of these quantities to their 
expected^ asymptotic behavior c n ~ (3~ n ~ l exp[B(n + I) 1 / 3 + O^n^ 1 ^ 3 )] with B = B x = 

3 ^ 2 /i 2 7,u/3 — in the case of the susceptibility, while B = B £ 2 = ^f^l in the case of 

^ 2 (/?). Both in the case of the sq and of the tr lattices, using the sets of coefficients of 
the susceptibility and of the correlation length, we can estimate 77 = 0.25(1) from the ratio 
B x /B^2. This estimate is however suspect because, assuming 77 = 0.25, we can consistently 
estimate b = 1.46(10) in the case of the sq lattice and b = 1.36(10) in the case of the tr 
lattice. Thus we must suppose that the series coefficients are not yet near enough to their 
asymptotic values that the value of b can be reliably estimated by this direct analysis, in 
spite of the fact that the ratio B x /B^ already takes the expected value. 

Finally, it is also interesting to notice that, as indicated in our previous analysis 1 , a study 
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of the function R(/3) = x(P) I £ 2_ri '(P) > °f its logarithm and of its log-derivative, supports 
the arguments of Ref. [9| concerning the absence of singular multiplicative corrections by 
powers of ln(£) to the critical behavior of x an d gives results completely consistent with the 
assumption that 77 = 1/4. Indeed if 77 7^ 1/4 or if 77 = 1/4, but there were in x singular 
multiplicative corrections, either Pade approximants or DAs of the above quantities should 
detect some singular behavior as (3 — > j3 c , which, at the present orders of expansion, is not 
seen at all. 

In conclusion, our recent HT series data confirm the BKT predictions for the critical 
behavior of the XY system and confirm or improve the existing estimates of the critical 
parameters for both the sq and the tr lattices. 
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FIG. 1: XY model on the sq lattice. Distribution of the singularities of a class of second-order 
DAs of ln 2 (x) vs their position on the (5 axis (open histogram). The central value of the open 
histogram is (3 C = 0.5599(7). The bin width is 0.0007. The vertical dashed line shows the critical 
value P c = 0.55995 indicated by the simulation of Ref. for which one can guess an uncertainty 
somewhat smaller than ours. The hatched histogram represents the distribution of the exponent 
a obtained from DAs of L{— c^X)^ biased with (3 C = 0.55995, vs their position on the a axis. The 
central value of the hatched histogram is a = 0.499(5) and the bin width is 0.003. 
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FIG. 2: A comparison of the estimates of the exponent a obtained from a class of DAs of ln 2 (x) 
and of L(— C2,x)> biased with the inverse critical temperature. Results are shown for both the sq 
and the tr lattices. We have varied the value of 0^ as , used to bias the DAs, in a small vicinity 
of the values (3l q = 0.55995, in the case of the sq lattice, and /3* r = 0.3412, in the case of the tr 
lattice. The temperature-biased exponent estimates are plotted vs x = 0^ as / ' @l q in the case of 
the sq lattice and vs x = [3^ as / '/3' r in the case of the tr lattice. The black squares (resp. black 
triangles) show the results obtained from the analysis of ln 2 (x) in the case of the sq lattice (resp. 

lattice) and the open squares (resp. open triangles) those obtained from the study of L(—C2,x) 
in the case of the sq lattice (resp. tr lattice). 
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FIG. 3: XY model on the sq lattice. A scatter plot of the singularities of a class of first-order 
DAs for L(—c S 2,x) m the complex (3 plane. Here x = Re((3) and y = Im((3). The central circle 
has radius (3 C . The small circles are drawn to enclose clusters of singularities which are likely to 
coalesce around tips of cuts. 
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FIG. 4: XY model on the tr lattice. Distribution of the singularities of a class of second-order inho- 
mogeneous DAs of ln 2 (x) versus their position on the [3 axis. The central value of the distribution 
is C = 0.3412(4). The bin width is 0.0003. 
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TABLE I: XY model on the tr lattice. The series expansion coefficients for the nearest-neighbor 
spin-spin correlation C(0, 0; 1, 0; (3), the reduced susceptibility x{P) an d the second moment of the 
correlation function m( 2 )(/3). 



order 


C(0,0;1,0;/J) 


x(P) 


m( 2 )(/3) 








1 





1 


1 


6 


6 


2 


2 


30 


72 


3 


7 
2 


135 


579 


4 


5 


570 


3834 


5 


35 
6 


2306 


22520 


6 


14 
3 


18083 
2 


121754 


7 


-81 
16 
-3769 
72 
-165161 
720 
-7821 
10 

-20160371 

OD^U 

4320 
-87289819 


276657 


4952033 


8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


8 

777805 
6 

14339641 
30 

208590287 

120 
8995595389 

Q1 QQ71 QQyc^ 

oiyy i too / o 

144 
65793037351 


8 

36001013 
12 

839474407 
60 

1264157753 
20 

400323755461 

1 AAf\ 
ouuo ( yiizoi i 
720 

12688548065393 


5040 
-10256893919 
226800 
-3357272555039 

29030400 
-1400375733941 

4838400 
-26431737035251 

37324800 
-55206137402197 
32659200 
-6827447251427903 


840 

165647319078571 

604800 
4600845479023849 

4838400 
1983863997387623 

604800 
24492996075345043 
2177280 
1671043059049640293 

43545600 
37817672635562705657 


2520 

3152231835174739 
151200 
411232110524237321 

4838400 
826713976281365323 

2419200 
4220638202244829129 

3110400 
662836750489256935 
124416 
2575741048252255298333 


1741824000 
-480824970393676609 


290304000 
4025832361031298767249 


124416000 
728769389306358221619671 


54867156000 


9144576000 


9144576000 
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